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ABSTRACT 

We have performed a pointed survey of N2D+ 2-1 and N2D+ 3-2 emission toward 64 N2H+-bright 
starless and protostellar cores in the Perseus molecular cloud using the Arizona Radio Observatory 
Submillimeter Telescope and Kitt Peak 12 m telescope. We find a mean deuterium fractionation in 
N2H+, Rd = A^(N2D+) / A^(N2H+), of 0.08, with a maximum Rd = 0.2. In detected sources, we find 
no significant difference in the deuterium fractionation between starless and protostellar cores, nor 
between cores in clustered or isolated environments. We compare the deuterium fraction in N2H+ 
with parameters linked to advanced core evolution. We only find significant correlations between the 
deuterium fraction and increased H2 column density, as well as with increased central core density, for 
all cores. Towards protostellar sources, we additionally find a significant anti-correlation between Rjj 
and bolometric temperature. We show that the Perseus cores are characterized by low CO depletion 
values relative to previous studies of star forming cores, similar to recent results in the Ophiuchus 
molecular cloud. We suggest that the low average CO depletion is the dominant mechanism that 
constrains the average deuterium fractionation in the Perseus cores to small values. While current 
equilibrium and dynamic chemical models are able to reproduce the range of deuterium fractionation 
values we find in Perseus, reproducing the scatter across the cores requires variation in parameters 
such as the ionization fraction or the ortho- to para-H2 ratio across the cloud, or a range in core 
evolution timescales. 

Subject headings: ISM: kinematics and dynamics - ISM: molecules - ISM: structure - radio lines: ISM 
- stars: formation 



1. INTRODUCTION 

Prior to the formation of a star, the cold (T < 20 K) 
and dense (n > 10^ cm~'^) conditions within star-forming 
molecular cloud cores drive a cold-gas chemistry that 
has been well-studied in recent years. Many molecular 
species, including CO and its isotopologues, become de- 
pleted in the gas phase in core centers by freezing out 
onto dust grains, forming an icy mantle. Depletion in 
CO has b een doc umented i n a numb er of source s (see , 
e.g. Kuipe r et al.llT9 96: Will acv et al. 1998; CascU i et al 



1999; .Bergin et al.ll2001 : Juvela et al.ii2002; .Bergin et al . 
2002 ). while iChristie et al. (2013) have performed the 



first systematic survey of CO depletion toward cores em- 
bedded in four nearby molecular clouds: Taurus, Ser- 
pens, Ophiuchus, and Orion. Being significantly absent 
from the gas, depleted species are thus no longer good 
tracers of core physical conditions, and the kinematics 
revealed by line profiles relate to the core outer envelope 
only. Molecular species which deplete at higher densities, 
such as the nitrogen-bearing species NH3 and N2H+, are 
better suited to probe the structure and kinematics of 
the interiors of dense cores. Many deuterated molecules, 
where one or more hydrogen atoms in a molecule is re- 
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placed by deuterium, are also preferentially formed in 
cold core interiors, making them excellent tracers of phys- 
ical conditions. 

In cold cores, abundances of deuterated molecular 
species can become enhanced relative to their undeuter- 
ated counterparts through reactions with deuterated H;^ , 
H2D+. H2D+ becomes mor e abundant at cold tempera- 
tures through the reaction ()Millar et al.lll989l ) 



H+ + HD o H2D+ + H2 + A£; 



(1) 



where AE = 232 K when all reactants are in their 
ground state. Equation [1] is thus exothermic in the for- 
ward direction, but AE is sufficiently small that the 
reaction proceeds in both directions at temperatures 
T > 20 K, resulting in no net increase in H2D+. In cold 
cores, however, the energy barrier is sufficient to signif- 
icantly lessen the backward reaction rate. H2D+ thus 
can quickly become abundant relative to H3", to levels 
orders of magnitude above the interstellar [D]/[H] ratio 
of - 10"^ (|York fc Rogersonl[l976h . The depletion of CO 
from the gas phase in cold cores removes an important 
H2D+ destruction pathway, allowing further enhance- 
ment of the deuterium fraction in Hj^. 

H2 and H2D+ exist in two states, ortho- and para-, 
depending on the spin of its two hydrogen atoms. The 
amount of H2D+ in the gas depends strongly on the 
ortho-to-para-H2 ratio. For example, the backwards re- 
action of Equation [T] proceeds more rapidly when the 
main reactants are ortho-H2D+ and ortho-H2. While 
the equivalent reaction with para-H2D+ and para-H2 
is endothermic and slow in cold regions, para-H2D+ 
is efficiently converted to ortho-H2D+ through reac- 
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tions with ortho-H2. The overall deuteration of Hg 
(and consequently of N2H+) will th us be reduced where 
the ortho-t o-para- H2 ratio is hi^h ([Pagani et al.l 119921 : 
iGerhch et al. 2002; iFlower et al.l l2006[ ). The deuterium 
chemistry is further complicated by additional parame- 
ters which may vary between star forming regions, in- 
cluding the grain size distribution (where d epletion onto 
gra ins decreases with increasing grain size; iCaselli et alJ 
|2008|) , and the ionization fraction in the dense gas (which 
leads to a greater dissociative recombination rate of 
H2D+). 

When H2D+ is abundant, additional reactions further 
propagate the deuteration abundance enhancement to 
molecules and molecular ions. Because of the need for 
high densities, cold temperatures, and depletion of CO, 
the enhanced deuterium fractionation of molecules typ- 
ically present in dense cloud cores is thus thought to 
indicate a dense core is more evolved, and perhaps more 
likely to be significantly centrally condensed and grav- 
itationally unstable. Indeed, chemical models predict 
the deuterium fractionation in N2H+, for example, is ex- 
pected to incre ase with time as the central density of the 
core increases (lAikawa et al.. 2005i). 

Observations of N2H+ and N2D+ toward star form- 
ing regions typically measure the deuterium fraction- 
ation, Ro, as the ratio of column densities, Rd ~ 
A^(N2D+)/iV(N2H+). Resuhs in low-mass starless aiid 
protostcllar cores ra nge between Rp < 0.1 (Tine et al 
2000; Tur neii 12001]: iCrapsi et al.l l200lj: iBeljoche et a! 

< 0.5 (iTine et~ 
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2006 
2000 

20071 : iBourke et al 



Caselli et al. 



2002; Cr apsi et al.ll2005t iDaniel et al 
I2OI2. ). Very high Rd ~ 0.7- 1 values 



have also been determined f or several objects, bot h in 
beam-averaged observations (iMiettinen et al.ll2012D. o r 
through radiative transfer modeling (jPagani et al. 

mm- 

Several surveys have looked for trends in the deuterium 
fractionation in N2H"'" with core parameters linked to 
evolutionary state, such as line width and line asym- 
metry, central density, temperature and level of CO de- 
pletion. In a stu dy of twelve low-mass pr otostcllar and 
five starless cores, [Roberts fc Miilml ((20071) found no sys- 
tematic difference in Ru between the two samples, and 
no significant trends in core physical parameters beyond 
a tende ncy for the st arless cores to have narrower line 
widths. ICrapsi et al.l ([2005) surveyed 31 low mass star- 
less cores, and found that while strong correlations were 
not found between Rjj and many parameters, there was 
a set of properties that, taken together, were associated 
with the most evolved cores. These 'evolved core' prop- 
erties include the column densities of N2H+ and N2D+, 
the deuterium fraction ation of N^ II^, and the amount of 
CO depletion. The Cr apsi et al.l sample included cores 
from several different nearby star forming regions, in- 
cluding the Perseus and Taurus molecular clouds. Anal- 
ysed separately, they found the subsample of Taurus 
cores presented a more homogenous set. In a study of 
twenty low mass prot ostcllar cores, selecte d frorn multi- 
ple molecular clouds, lEmprechtinger et al.l ([2009) found 
correlations between the deuterium fractionation with 
dust temperature and CO depletion, where colder, more 
CO-depleted cores had greater Rd values. The authors 
also showed that cores with high Rd values were likely 
to show infall asymmetries in their spectra. These trends 



were strongest, however, when the sample was limited to 
the thirteen sources in the Perseus molecu lar cloud. In a 
sample of dark clouds. [Daniel et al.l ()2007[ ) show that the 
deuterium fractionation in N2H+ is efficient at densities 
n > 2 — 3 X 10^ cm~^ . In contrast to the tre nds found by 
lEmprechtinger et al.l for the Perseus cores, IDaniel et al.l 
found no trend in the column density ratio of N2D"'" 
to N2H+ with temperature in their homologous sam- 
ple. Comparing results from studies focusing on cores in 
single regions versus those in different environments sug- 
gests that the external environment plays a role in the 
chemical evolution of star-forming cores. 

Here, we present the first systematic pointed study 
of the deuterium fractionation in N2H+ toward a large 
sample of starless and protostcllar cores within a single 
molecular cloud, Perseus, to characterize those factors 
which impact the core chemistry while excluding cloud- 
to-cloud variation. We have observed N2D+ 2-1 and 
3-2 toward 64 dense cores in the nearby (250 pc; e.g., 
lCernisllT993[ ) Perseus molecular cloud, selecte d based on 
their N2II+ 1-0 brightness in the catalogue of lKirk et al.l 
(|2007| ). In addition, many of the targets have been pre- 
viously s urveyed in infrar ed (Evans ct al. 200tt). submil- 
limeter (|Kirk et al.l 120061 ) and millimeter ([Enoch et al.l 
[2OO6I) continuum emission , as well as in multiple molec- 
ular lines, inclu ding CO (iKirk et all 120071 : iRidge et al] 
[2OO6I) and NH3 (IRosolowskv et al.ll200 8D. 

We present the source selection and observations in fj2l 
and describe the spectra and analysis in [JH While we 
find several significant correlations between the N2H+ 
deuterium fraction with parameters expected to be indi- 
cators of evolved star-forming cores, including with II2 
column density and core central density, there is signif- 
icant scatter in the trends. Other parameters, such as 
the kinetic gas temperature, Tx, show no correlation 
with Rd, in contrast with the results of several previ- 
ous studies. Some of the observed scatter may be due 
to variations in the filling factor of the N2D+ emission 
within our beam, particularly for protostcllar sources, 
but we show that this effect should be relatively small. 
On average, we find that the Perseus cores have low CO 
depletion, but the few cores with high CO depletion also 
show higher N2H+ deuterium fractions, as expected from 
chemical models. For the protostcllar cores, we see a 
significant anti correlation between the bolometric tem- 
perature of the protostar and the N2II+ deuterium frac- 
tion. In addition, the protostcllar sources least likely to 
be detected in N2D+ emission also have low envelope 
masses, showing that the deuterium fraction in proto- 
stcllar sources decreases as protostars age, heating their 
surrounding envelopes, and losing envelope mass through 
radiation, stellar winds, and molecular outfiows. In [|4|we 
discuss the deuterium fractionation results in the Perseus 
cores with respect to the physical properties of the cores 
and embedded protostars, and compare our findings with 
chemical models, finding that current equilibrium and 
dynamic chemical models are able to reproduce the range 
of deuterium fractionation values we find in Perseus, but 
require variation in parameters such as the ionization 
fraction or the ortho- to para-H2 ratio across the cloud, 
or a range in core evolution timescales, to reproduce the 
scatter across the Perseus cores. We summarize our re- 
sults in iJSj 
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Y 


150 


3 


25 


25 


7 


30 


45 


01 


6 


4.05 





49 


9 


4 


5 


4 


(0 


0) 


85 


0.83 




152 


3 


25 


22 


5 


30 


45 


06 


5 


4.08 





50 


10 


9 


5 


8 


(0 


2) 


86 


1.64 


Y 



Note. — Units of right ascension are hourSj minutes, and seconds, and units of declination are degrees, arcminutes, 
and arcseconds. All core ID num bers and N2H"*" em ission line properties are taken from Kirk et al. (2007). The SCUBA 
data were originally described in lKirk et al.l II2006I ). and have been convolved to a 30" FWH M for this analysis. The 
designation of a core as protostellar relies on multiple literature sources, described in detail in ^2.11 

^ N2H"'" core ID numbers. 

^ N2H+ fjsri FWHM Av, T X (Tei; — Tjg), and T^^ based on hyperfine spectral line fitting. The tdv function was used 
to calculate r X {T^^ ~'^bg)- 

" SCUBA ID number and total SCUBA flux within a 30" FWHM beam. 

Protostellar designation as discussed in the text. 
° Sources where two velocity components along the line of sight are identified in N2H+ 1-0 line profiles. 
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2. SOURCE DATA 
2.1. Source selection 

The observed N2D+ cores wer e selected based on their 
bright N2H+ 1-0 emission from iKirk eFall ()2007[ ). but 
nevertheless span a wide rang e of condi tions in the 
Perseus molecular cloud. The IKirk et all sample in- 
cluded candidate dense cores from submillimeter data, 
points of high visual extinction from Palomar plates, 
and peaks of large-scale extinction from 2MASS data. 
Within our subset of 64 sources, fifty-one are located at 
the peak of an 850 /im continuum clump. Some core 
properties (Right Ascension, Declination, velocity of the 
local standard of rest, vlsr, FWHM line width, Av, and 
integrated intensity of associated N2H+ 1-0 emission, as 
well as SCUBA fiux and identification as protostellar or 
starless) are listed in Table [1] Based on submillime- 
ter continuum measurements, the cores in our sample 
have masses ranging from ~ O.IM0 to ~ lOM© (me- 
dian AI — 1.4 A/0), effective angular radii from ~ 15" to 
~ 55" (median r = 28"; 0.03 pc or 7000 AU at 250 pc), 
and have a medi an average density of ^ 2 x 10^ cm~'^ 
(|Kirk et al.ll20M) . 

To identify cores as starle ss or proto s tellar, we be- 
gin with the classification of IKirk et all (|2007D , which 
is based on analysis of Spitzer data in Perseus by 
iJorgensen et al.l ()2007| ). We further compare our target 
list with the catalogue of lEnoch et al.l ()2009t ). who iden- 
tified protostellar sources using Spitzer and 1 mm Bolo- 
cam data, and classify as protostellar those cores which 
have a protostellar source within 20" of the pointing lo- 
cation. Three additional Perseus cores in our sample 
show evidence of harboring embedded protostellar ob- 
jects in recent studies based on molecular line o r deep 
continuum observa tions (Schncc et al. 2012; Che n et all 
l2010tlEnoch et al.|[20iaiSchnee et al.. ,2010'). despite hav- 
ing been previously classified as 'starless' in the studies 
described above. Including these sources in our pro- 
tostellar subset, we find that a greater fraction of the 
Perseus cores are classified as protostellar (43/64) than 
starless (21/64). We further classify cores as 'clustered' 
if they reside within the NGC 1333 or IC348 regions, 
foll owing the boundary de finitions of each region given 
bv |j0rgensen et all (|2006t ). The targets are evenly di- 
vided between those found in the clustered star-forming 
regions (32 / 64) , and those in more isolated parts of the 
cloud (32 / 64) . The more isolated cores nonetheless reside 
largely in well-known groupings within Perseus, includ- 
ing the B5, Bl, and L1448 regions. 

2.2. SMT observations of N2D+ 3-2 

Sixty-four Perseus cores were observed at the Arizona 
Radio Observatory Submillimeter Telescope (SMT) in 
2010 April 10-11, November 27-28, and 2011 February 16. 
During the April epoch, the observations were performed 
using the ALMA Band 6 prototype receiver (1 mm) with 
the Chirp Transform Spectrometer (CTS) as the back- 
end. The spectral resolution was 47 kHz, corresponding 
to 0.06 km s~^. Most sources were observed for three 4 
minute scans in position-switching mode, with additional 
integration time on several weak sources to increase the 
signal to noise ratio. Pointing observations were per- 
formed every 1-1.5 hours by offsetting the detector fre- 
quency to the i^CO 2-1 transition (230.5 GHz). Dur- 



ing the later two sets of observations, the CTS was not 
available, and the filterbank backend was used instead. 
The spectral resolution was 250 kHz, corresponding to 
0.32km s~^. All sources observed with the spectrome- 
ter were re-observed with the filterbanks. Each source 
was typically observed for two 5 minute scans on both 
the November and February runs, with several final ad- 
ditional scans on several low signal-to- noise (S/N) detec- 
tions in the November run. The SMT beam at 231 GHz 
is -30" (FWHM), well-matched to both the N2H+ (26" 
FWHM) and NH3 (32" FWHM) data. 

All spectra were reduced in CLAS^. We first flattened 
each spectrum using a 4th order polynomial. Several 
scans had highly irregular baselines which could not be fit 
with the polynomial; these were discarded from further 
use, and, where possible, re-observed. 

The spectra were calibrated for each observing run 
based on observations of Mars (data taken in April) 
and Jupiter (data taken in November and February). In 
November, the vertical (V) polarization data showed in- 
consistent results in calibration scans, and these data 
were omitted from the final dataset. The calibration 
efficiencies were found to be 0.82 ± 0.07 for the April 
run, 0.68 ±0.02 for the horizontal (H) polarization of the 
November run, and 0.55 ± 0.02 and 0.48 ± 0.02 for the 
H and V polarizations of the February run. The change 
in calibration efficiency between November and Febru- 
ary may have been caused by a warmup of the receiver 
between the two epochs. Several Perseus targets with 
strong N2D"'" detections were observed in more than one 
run, to allow for a secondary check of the beam efficien- 
cies derived. The relative efficiencies were found to be in 
agreement with the planet-derived values, within uncer- 
tainties of '~ 15 % in the observed peak line intensities. 

After applying the calibration factors, all spectra of 
each source were summed together. We show in Figure 
[1] a comparison of the N2D+ 3-2 spectra toward four 
sources which were observed with both the CTS and fil- 
terbank backends. The Figure shows that the filterbank 
data recover the line emission for both faint and strong 
emission lines, but sources with narrow lines will suffer 
from artificial line broadening and depressed peak line 
brightnesses due to the poor spectral resolution of the 
filterbanks. In particular, cores with both faint and nar- 
row N2D+ emission that would be detected by the CTS 
observations may not be detected in the filterbank data 
(see Pers 3 in the Figure). Of the eleven sources observed 
with both the CTS and the filterbanks, only one source 
(Pers 3) showed line emission in the CTS data but not in 
the filterbank data. For this source, we include only the 
CTS observations in the final data. Our analysis below 
depends only on the N2D+ 3-2 line centroids and inte- 
grated intensity, which agree well between the CTS and 
filterbank data, and we thus conclude that our results are 
not significantly affected by the low spectral resolution 
of the data. 

The N2D+ 3-2 emission line contains hyperfine struc- 
ture, but no structure was visible at a 3(7 level above the 
rms noise in the data for any Perseus sources (Pers 25 
shows a marginal detection of the satellite lines). A sin- 
gle Gaussian fit in CLASS was thus used to characterize 

^ CL ASS is part of the GILDAS software suite available from 
|http://ww w.iram.fr/IRAMFR/GILDAS 
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Fig. 1. — Comparison of N2D"'" 3-2 spectra obtained toward four Perseus sources with the CTS (black hnes; Kght grey Knes show the 
CTS data smoothed to the filter bank resolution) and the filterbank backends (grey lines; red in the online version). 



the emission, and the returned area from the fit is re- 
ported as the line integrated intensity, W32, for sources 
detected with S/N > 3. To test the effects on our re- 
suhs from the low velocity resolution of the filterbank 
data, and the fitting of a single gaussian rather than the 
full hyperfine structure to the lines, we additionally per- 
formed hyperfine fits toward those sources observed at 
high velocity resolution with the CTS. We find that the 
line centroids are recovered by the single gaussian within 
the filterbank velocity resolution of 0.32 km s~^, but the 
returned gaussian line widths are a factor ~ 2 larger than 
those found by the hyperfine fits. The returned area from 
the gaussian and hyperfine fits additionally agree within 
uncertainties. Final sensitivities for each source are 
listed in Table H 

2.3. Kitt Peak 12m Telescope observations of N2D'^ 

2-1 

Observations of N2D+ 2-1 were made toward 54 
Perseus cores within the allotted time at the Arizona 
Radio Observatory Kitt Peak 12 m over 2010 Novem- 
ber 29 and 30, and December 11, 12 and 13. Using the 
2 mm receiver, the observation rest frequency was set to 
151.217178 GHz, the frequency of the strongest hyperfine 
component of the 2-1 transition. The backend used was 
the Millimeter Autocorrelator, with 100 MHz total band- 
width and 6.1kHz spectral resolution (0.012 km s~^ at 
151.2 GHz). Sources were observed in position-switching 
mode for two or three periods of six minutes each to 
achieve similar rms noise levels based on weather condi- 
tions. The data were corrected at the telescope for the 



atmosphere and for telescope losses. The aperture effi- 
ciency of the telescope is 45% at 145 GHfl Further cali- 
bration to the Tmb scale was done using observations of 
Jupiter, which were performed once each shift. Based on 
the Jupiter measurements, uncertainty in the final cali- 
bration scale is ^ 5%. At 151 GHz, the telescope beam 
is - 41" (FWHM). 

Spectra were reduced in CLASS. The data were 
smoothed to a final velocity resolution of 0.14kms~"'^. A 
linear baseline was removed and the data were summed 
for each source. Detections were fit using the hyperfine 
structure (hfs) CLASS routine, usin g line frequen c ies an d 
relative strengths calculated from iPagani et al.l ()2009D . 
The resulting line centroid velocities, vlsRi and widths. 
At;, should have better accuracy than would be provided 
by a single Gaussian fit that did not take into account 
the hyperfine structure of the line. We do not discuss fur- 
ther the derived N2D"'" 2-1 centroid velocities, however, 
due to an unexplained, episodic velocity shift observed 
in some of the data. While this observed velocity shift in 
the observed line, when present, was of a similar magni- 
tude in nearly all clearly detected sources, there was no 
clear pattern in time when the offset occurred, and no 
obvious cause in the telescope log. We find no evidence 
in N2D"''-bright sources that the velocity shift occurred 
between scans of the same source, so are confident the 
resulting summed data accurately represent the source 
line emission. In nearly all sources, the returned line 
opacities also have large uncertainties due to a lack of 

* http:/ /aro. as. arizona.edu/12_obs_manual/chapter_3. htm 
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high S/N in the satelhte components. We therefore pri- 
marily use the hfs-fit vefocity information to determine 
three windows over which we measured the integrated 
N2D+ 2-1 intensity for each core, W21, based on the hne 
vlsr and Av, and the expected focations of N2D+ 2-1 
hyperfine components. Where no hne was detected, an 
upper hmit to the integrated intensity was determined, 
foUowing the same method but defining spectral windows 
based on the centroid velocity and width of the observed 
N2H+ 1-0 line toward the core. Final sensitivities for 
each source are listed in Table [2] 

2.4. Other data 

In our analysis, we make use of several published cata- 
logs focused on the Perseus molecular cloud. The N2H+ 
1-0 catalog used to select cores for N2D+ observations 
was based in part on previous 850 /im continuum obser- 
vations performed with the James Clerk Maxwell Tele- 
scope using the Submillimetre Com mon User Bolome - 
tcr Array (SCUBA) at 15" FWHM (|Kirk et al.ll20M ). 
smoothed to a final effective resolution of 20" FWHM 
to reduce pixel-t o- pixel noise. In their N2H+ study, 
iKirk et al.l 1)200 7') also present observations of C^^O 2-1 
made with the IRAM 30 m telescope at 11" FWHM. AU 
cores with N2H+ emission were also detected in C^'^^O. 
The C^*0 emission was comprised of multiple veloc- 
ity compon ents for many cores, and in their analysis, 
IKirk et all determine the individual C^*0 component 
that best matches the N2H+ emission for each core. We 
use these matched components in our analysis and dis- 
c ussion below. 

iRosolowskv et all (|200 8l performed a targeted study of 
NH3 (1,1) and (2,2) emission toward a large number of 
sources in Perseus, calculating kinetic temperatures for 
each core. Their targets were selected from several cat- 
alogs, foremost a 1.1mm study of the region, observed 
using BOLOCAM at the Caltech Submillimeter Observa- 
tory (Enoch et al. 2006). Consequently, the NH3 obser- 
vation centers are frequently offset relative to our N2H+ 
and N2D+ observations. For our targets, however, most 
of the offsets remain small (< 15", or half the N2D+ 
3-2 beam FWHM at the SMT). Of those Perseus cores 
detected in N2D+ emission with offsets of ^ 20"- 25" 
between the NH3 and N2H+ targets (five cores), only 
one shows a moderate (~ 0.5kms~^) offset in line-of- 
sight velocity between the NH3 and N2H+ observations, 
with the rest in agreement within observational uncer- 
tainties. We thus expect that the N2H+ and NH3 ob- 
servations are tracing similar material for most of our 
targets, and that the kineti c temperature me a surem ents 
from the NH3 analysis bv IRosolowskv et"aLl ()2008D are 
accurate for the cores observed in our study. Four cores 
in our catalogue (23, 71, 125, and 148), which were all 
detected in N2D+ emission, have no nearby counterparts 
in the NH3 dataset. 

Fig. Set 2. Spectra of N2H+ and N2D+ 

3. RESULTS AND ANALYSIS 

We list the parameters of the fits to detected sources 
in Tabled] For this study, a core is 'detected' if emission 
was seen with S/N > 3 in Tmb in either the N2D+ 2-1 
or 3-2 transitions. Of the 54 Perseus cores observed in 
emission from both N2D+ 2-1 and 3-2, 45 were detected 






J\J.^//W\rU\^rr/yWJ^•/''-'^■'^^^ 



km s 

Fig. 2.— Spectra of N2H+ and N2D+ toward each of the 64 
Perseus cores in this study. N2H+ 1-0 spectra are shown in black, 
N2D+ 2-1 spectra are shown in grey (blue in the online version), 
and N2D+ 3-2 spectra are shown in dark grey (red in the onUne 
version). For clarity, the N2D+ 3-2 data have been multiplied by a 
factor of 2, and spectra are offset from zero to display all available 
datasets. Some cores were observed in one N2D+ transition only. 
The N2H+ 1-0 data were first presented in Kirk et al. (2007). 



in both transitions with S/N > 3. Two sources were 
detected only in N2D+ 2-1 emission, four were detected 
only in N2D~'' 3-2 emission (likely due to insufficient sen- 
sitivity in the N2D+ 2-1 observations), and three were 
undetected in either transition. Of the remaining ten 
cores which were not observed in N2D"'" 2-1 emission, 
nine were undetected in emission from N2D'^ 3-2. We 
show in Figure the N2H+ 1-0, N2D+ 2-1, and N2D+ 
3-2 spectra toward each target. 

We can further break down the detections and non- 
detections as a function of whether the core has pre- 
viously been identified as starless or protostellar, and 
whether the core lies within a clustered region of the 
Perseus molecular cloud. Based on these definitions, 36 
of 43 (84 %) protostellar sources were detected, and 16 
of 21 (76 %) starless sources were detected. Of the clus- 
tered sources, 29 of 32 (91 %) were detected, and 23 of 
32 (72 %) isolated sources detected. 

3.1. Kinematics 

The centroid velocity of a core is an easily measured 
property. As described above, even with relatively low 
velocity resolution, we are able to determine the N2D+ 
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TABLE 2 

N2D+ LINE PROPERTIES 



ID 



■VLSR 
km 



Ad 
km 



N2D+ 
Tr 
K 



2-1 



rms 
K 



W21_ 
K km s ^ 



VLSR 
km 



N2D+ 3-2 

At) Ipeak 

kms~^ K 



rms 
K 



Ws2 
K km 



10.38 
8.91 
8.36 
9.06 
8.60 
7.78 

10.05 
7.28 
6.29 
7.04 
6.89 
7.05 
6.67 
7.55 
6.98 
6.10 
8.79 
8.13 
8.12 
7.68 
8.09 
7.76 
7.72 
9.20 
7.68 
8.27 
8.48 
8.49 
6.90 
8.28 
7.90 
7.91 
8.20 
7.29 
7.88 
8.56 
8.71 
7.66 
4.86 
3.91 
6.02 
3.61 
4.39 
3.52 
3.09 
3.15 



0.04 
0.05 
0.03 
0.03 
0.02 
0.06 
0.03 
0.04 
0.04 
0.04 
0.07 
0.04 
0.04 
0.03 
0.04 
0.02 
0.06 
0.02 
0.11 
0.04 
0.03 
0.03 
0.02 
0.05 
0.02 
0.03 
0.02 
0.05 
0.04 
0.03 
0.05 
0.07 
0.06 
0.04 
0.04 
0.04 
0.03 
0.05 
0.06 
0.48 
0.05 
0.03 
0.02 
0.02 
0.03 



0.46 
0.34 
0.87 
0.39 
0.16 
0.34 
0.94 
0.94 
0.63 
0.38 
0.77 
0.24 
0.77 
0.30 
0.31 
0.36 
0.63 
0.36 
0.68 
0.31 
0.66 
0.85 
0.37 
1.28 
0.99 
0.35 
1.09 
0.93 
0.44 
0.64 
0.47 
1.35 
0.75 
0.40 
0.56 
0.53 
0.41 
0.52 
0.49 
2.05 
0.51 
0.73 
0.73 
0.46 
0.51 



0.09) 
0.09) 
0.01) 
0.07) 
0.06) 
0.12) 
0.08) 
0.13) 
0.11) 
0.07) 
0.16) 
0.08) 
0.12) 
0.06) 
0.08) 
0.05) 
0.10) 
0.04) 
0.25) 
0.09) 
0.07) 
0.13) 
0.04) 
0.20) 
0.05) 
0.06) 
0.07) 
0.22) 
0.09) 
0.09) 
0.14) 
0.16) 
0.20) 
0.09) 
0.10) 
0.15) 
0.08) 
0.09) 
0.10) 
0.74) 
0.10) 
0.10) 
0.05) 
0.05) 
0.06) 



Detected in both N2D+ 2-1 and 3-2 
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line centroid with good accuracy. Previous work has 
shown that the velocity centroids measured for denser 
core gas, as traced by, e.g., N2H+, tend to be similar to 
those measured for the surrounding, lower density enve- 
lope of ma t erial, as traced b y, e.g., C^®0 (jWalsh et al.l 
[2OOI [20071 [Kirk et al.l 120071) . For cores in the Perseus 
molecular cloud, the differences in centroid velocities be- 
tween N2H+ and C^^O are small, typically less than the 
sound speed, but are often non-zero (rms absolute veloc- 



ity difference 0.17 km s^ MKirk et al.ll2007D . A much 
tighter relationship has been found for the centroid ve- 
locities of N2H^ and NH3 in the Perseus cores, whe re the 
NH3 data was taken from R osolowskv et al.l ()2008D . The 
mean absolute offset of N2H"'" to NH3 centroid velocities 
is 0.07kms~^ (jJohnstone et al.ll2010f ). This measure- 
ment suggests that the N2H+ and NH3 occupy a similar 
volume within the dense core. 
Here, we compare the centroid velocities measured us- 
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Fig. 3. — Left: The distribution of centroid velocity differences for N2H+ to NH3 and N2H+ to C^^O. The vertical dotted lines denote 
the approximate sound speed in the region (assuming T=15 K), and km s~^. Right: The distribution of centroid velocity differences for 
N2D+ to N2H+, N2D+ to NH3, and N2D+ to C^*0. The vertical dotted lines denote the approximate sound speed in the region (assuming 
T=15 K), and km s"!. 



ing the N2D+ 3-2 transition to the various other hne cen- 
troids aheady measured for Perseus dense cores. Since 
N2D+ is expected to trace a higher density of gas than 
even N2H+, it could be expected to show an equal or 
larger typical velocity difference to the lower density gas 
tracers than N2H+. The left panel of Figure [3] summa- 
rizes the previously observed velocity differences between 
N2H"'", NH 3 and CO in the Perseus cores described in 
iKirk et aD (|2007j l. here including only cores which are 
part of this survey. The dashed blue line shows the 
velocity di f ferenc e in N2H+ and C^^O observed in the 
IKirk et al] ()2007D IRAM pointed survey, while the pur- 
ple line shows the velocity difference between N2H+ and 
NH3 . The right panel of Figure [3] shows the comparable 
velocity differences found using N2D+. Similar to the 
trend seen in N2H+, much larger centroid velocity differ- 
ences are seen between N2D+ and C^^O than N2D+ and 
NH3. The velocity differences between N2D+ and N2H+ 
appear to be slightly larger than for N2D+ to NH3, al- 
though statistical tests (see below) show that there is no 
significant difference between them. Although the C-'^^O 
data were taken with significantly better angular resolu- 
tion than the N2H+ and N2D+ observations, the emis- 
sion traces the outer regions of the core due to the effects 
of depletion and the lower critical density of the C^^O 
transition. C^*0 observations made with a larger beam 
would include more emission from lower density material, 
and would likely show larger velocity differences between 
the tracers. 

We next compare the velocity difference distributions 
for N2H"'" versus NH3 and C^^O for all of the dense cores 
detected in N2H+ versus the subset observed and de- 
tected in N2D"'" , using a two- sid ed Kolmogorov-Smirnov 
(KS2) test (eg, iConoveil ITqM ). The KS2 test shows 
the likelihood that two distributions are drawn from the 
same parent population. The KS2 probabilities were 
found to be roughly 90% or above, indicating that the 
cores m which N2D+ was detected are representative, in 
terms of centroid velocity differences, of the full N2H+ 
sample of cores. Furthermore, we compare the distribu- 
tion of velocity differences seen between N2D+ or N2H+ 
and NH3 or C^^O, and found that the KS2 probabilities 
are sufficiently large (roughly 50% or above) to suggest 
that N2D+ and N2H^ behave similarly with respect to 



other gas tracers. 

We also search for trends in the velocity difference be- 
tween N2D+ and N2H+ in the protostellar and starless 
core populations, as well as in the clustered and isolated 
cores. A comparison between the velocity differences 
for starless and protostellar cores shows that the former 
has a slightly smaller dispersion [(T(Ay) — 0.07f kms^^ 
versus 0.09f kms^^]; however, a KS2 test suggests that 
these differences are not statistically significant. Simi- 
larly, splitting the cores into those which are found in 
more isolated conditions and those found in more clus- 
tered environments shows that the former tend to have 
smaller velocity differences [cr(Ay) = 0.067 km s^^ ver- 
sus 0.097 km s""'^], but again, a KS2 test shows these dif- 
ferences are not statistically significant. 

3.2. N2D^ column density and the N2H^ deuterium 
fractionation in Perseus cores 

Both the N2D+ 2-f and 3-2 transitions contain hyper- 
fine structure, which allows independent calculation of 
the excitation temperature. Tea;, and the opacity, r, of 
the emitting gas. In practice, however, a high signal- 
to-noise ratio is needed in the satellite components to 
calculate T^x and r with precision. No N2D"'" 3-2 spec- 
tra showed significant hyperfine structure above the rms 
noise level. While hyperfine structure in the N2D"'" 2-f 
line is seen toward multiple sources, our hyperfine line 
fits are unable to constrain well the N2D"'" 2-f line opac- 
ity and the resulting uncertainties on T^x and r are large 
(see Tabled]). In contrast, the observed cores were se- 
lected t o be brigh t in N 2H"'" f-0 emission, and conse- 
quently EkkeTd] (|2007[ ) were able to fit well the N2H+ 
hyperfine structure. The N2fl''" excitation temperature 
can be determined from their fit results. 

We expect that N2H+ and N2D+ trace similar phys- 
ical conditions given their similar kinematics in the 
Perseus cores (see i)3.ip . While we are unable to di- 
rectly compare the excitation of the lines in our cores, 
previous studies have s hown that the exc i tation of the 
two s pecies are similar (jGerin et al.l l200fl : iCaselli et al.l 
l2002f ). We can then test our assumption that the N2D+ 
lines are optically thin by performing a hyperfine fit to 
the line structure, setting T^x equal to that derived from 
the N2H+ f-0 analysis. For both N2D+ transitions, the 
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total optical depths thus determined are universally < 1. 
We are thus confident that we can assume the N2D"'" 
emission is optically thin in the following analysis. 

In the optically thin regime, we calculate the N2D+ 
column density, A^(N2D+), from the N2D+ integrated 
intensity W following the form ula for optically thin emis- 
sion from lCaseUi et al.l (|2002f ): 



iV(N2D 



_ SttW gi 



1 



gu Jv{Tex) - Ji^{Tbg) 

1 Qrot 

1 - exp{-hiy/kTex) giexp{-Ei/kTe, 



(2) 



Here, A is the wavelength of the observed transition, 
A is the Einstein coefficient {A — 1.97 x 10^"* and 
7.14 X lO"'^ for N2D+ 2-1 and 3-2, respectively), gi 
and gu are the statistical weight of the lower and up- 
per levels, h is Planck constant, Tex is the excitation 
temperature, Qrot is the partition function, Ei is the 
energy of the lower level, and k is the Boltzmann con- 
stant. Jv{Tex) and J,y{Thg) are the equivalent Rayleigh- 
Jeans excitation and background temperatures. For 
linear molecules like N2D+, Ei = J; (J; + l)hB and 
Qrot = YALf)gi^'y^v{-Ei/{kTex)), where J; is the rota- 
tional quantum number of the lower level, and the rota- 
tional constant B = 38554 x 10^ Hz. For each core, we 
use the N2D+ 3-2 integrated intensity and set the exci- 
tation temperature, T^x to that derived from the N2H+ 
1-0 analysis. For sources where N2D+ 3-2 was not de- 
tected, we instead calculate the N2D+ column density 
from the N2D+ 2-1 integrated intensity. Where neither 
N2D+ line was detected, we determine an upper limit 
to the integrated intensity based on the rms noise of the 

1 /2 

spectra, aw, following aw = H^^VresCTrmsN^l^^^. Here, 
Awres is the spectral resolution, arms is the 1—a noise 
value, and Nchan is the number of channels over which 
the line integrated intensity was determined. The N2D+ 
column density upper limit is then calculated using the 
N2H+ 1-0 Tex for the source and setting = 3 x aw 
for each of the N2D+ 2-1 and N2D+ 3-2 transitions, and 
selecting the resulting lowest limit. 

The deuterium fraction in a species, Rd, is defined 
as the ratio of the column density of the deuterated 
molecule with that of its hydrogen-bearing counterpart, 
where Rd = A^(N2D+) / A^(N2H+) for N2H+. For 
sources which show no evidence of multiple velocity com- 
ponents along the line of sig ht, we calculate A^(N2H+) 
following iCaselli et al.l ()2002D for optically t hick transi- 
tions, using the hyperfine fit parameters in IKirk et all 
(|2007f) . Two velocity components are visible in the 
N2H"'" spectra for nine cores, but we are unable to sep- 
arate velocity components in the N2D+ spectra. To 
provide a more accurate estimate of the N2H"'' deu- 
terium fraction for these objects, we thus calculate the 
N2H+ column density using the tot al integra ted inten- 
sity of the N2H+ line pre sented in IKirk et al. (2007), 
following Equation A12 in iJohnstone et al.l ~f2010') and 
assuming optically thin e mission. Where targets over- 
lap with IJohnstone et al.l . our N2H+ column densities 
agree within uncertainties. The uncertainty in Rd is 
propagated from the uncerta inty in both iV(N2D^) and 
Af(N2H+) (taken to be 20 %; 'Johns tone et aLllMO) . 

The resulting deuterium fractionation of N2H+ in each 



core is listed in Table [31 We find a mean Rd = 0.08 for 
all cores, with a standard deviation of 0.04. We plot 
7V(N2D+) against iV(N2H+) in Figure g] (left), overlaid 
with lines of constant Rd — 0.01, 0.08, and 0.25. Cores in 
which one or both N2D+ lines were detected are shown 
separately, and upper limits are shown for those cores 
where neither line was detected. Single-line detections 
are more likely at lower N2D"'" column densities, but span 
a range in N2H+ column density, and hence a range in 
Rd- 

In Figure [4] (right), we plot Rd &s a, function of 
A^(N2H+). No separation is seen in Rd between proto- 
stellar and starless sources. Intriguingly, the highest Rd 
values are found toward cores with relatively low N2H"'" 
column densities compare d with the ent i re sam ple. A 
similar result was found by IFontani et al.l (|2011| ) in their 
N2H+ / N2D+ survey of massive star forming regions, 
which they suggest may be due to a decrease in the deu- 
terium fractionation, while the N2H+ abundance contin- 
ues to increase, in the most embedded protostellar phase. 
Figure |4] (right) shows that on average, the protostellar 
targets in Perseus do have greater N2H+ column densi- 
ties than the starless cores. The high Rd, low A^(N2H+) 
points include both starless and protostellar cores. We 
will show in [21 however, that these protostellar sources 
are consistent with being young. 

We find that there is little difference in the deuterium 
fraction Rd values between the starless and protostellar 
core subsamples, nor between the isolated and clustered 
core subsamples. All subsets have similar mean values 
and standard deviations to that of the overall sample. 
This can be seen clearly in Figure |S1 where we show his- 
tograms of the Rd values for the Perseus cores, separated 
by subsample. Most of the Perseus cores have deuterium 
fractionation values Rd ^ 0.1; however a small number 
of cores form a tail in the distribution to higher Rd val- 
ues. As the Figure shows, cores with higher deuterium 
fractionation are not confined to the starless or proto- 
stellar population, nor are they found only in clustered 
or distributed star forming environments. 

To further investigate the properties of the high-i?/) 
cores, we show in Figure |6l the distribution of Rd values 
across the Perseus molecular cloud, colored to show cores 
with Rd > 0.1, Rd < 0.1, and N2D+ non-detections. 
Most of the well-known groupings of star forming cores 
(e.g., B5, IC348, NGC 1333 and L1448) contain at least 
one core with Rd > 0.1, as well as a few moderately 
deuterated cores and N2D"'" non-detections, and we see 
no systematic trend of the N2H"'" deuterium fraction 
across the cloud as a whole. Five cores have Rd > 0.15. 
Comparing the physical properties of these cores with the 
rest of the sample, we find that the high Rd cores are 
characterized by greater 850 /im continuum flux and core 
concentration, larger column densities of H2, as well as 
higher NH3 and lower N2H+ column densities (as seen 
in Figure HI (right) . 

Similar results have been found previousl y in Perseus. 
For example, toward three cores in Perseus, iRoueff et al.l 
(2005) detected triply-deuterated ammonia, ND3, with 
an abundance of 0.1% relative to NH3, suggesting the 
abundance of singly-deuterated NH3 is ~ 10%. The 
average Rd values found here also agree w ell with 
other studies of low-mass cores. iCrap si et al.l find a 
mean Rd = 0.1 for their core sample, with values 
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Fig. 4. — Left: Af(N2D+) vs. Af(N2H+) for the Perseus cores. Cores with both N2D+ 3-2 and N2D+ 2-1 detections are shown by solid 
black triangles, while cores with a single detection in either N2D+ 3-2 or N2D+ 2-1 are plotted with open black triangles. Arrows show 
the Af(N2D"'") upper limits for cores with non-detections in both transitions. The dotted lines represent constant N2H+ deuterium fraction 
levels, Rd = 0.01, the mean Rd = 0.08 for the sample, and Rd = 0.25. Right: Ro vs. 7V(7VN2H+) for the Perseus cores. Starless and 
protostellar cores are shown by black and grey (colored red in the online version) triangles, respectively, while deeply embedded sources 
are shown as open triangles (colored blue in the online version). 



ranging from 0.03 to 0.44. Similarly, other studies of 
low mass starless and protostellar cores find Rd val- 
ues wit hin the range of a few hundredths to a few 
tenths (iDaniel et al.l IgqOTt lEmp rcchtingcr c iraLl 120091: 
Roberts fc Millaij I2007t iFriesen et al.ji2010t iGerin et al l 



2001; Micttine iTet all 12009: Pa gani et an " l2009[ ). Iiiter 

estingiy, no core in Perseus shows N2H+ deuterium frac- 
tionation levels above ~ 0.25, despite our core sample 
spanning mean and central densities up to n ^ 10^ cm~^. 
Regions of high deuteration exist in Perseus, but are 
not apparent in our a n alysis . In NGC 1333, for exam- 
ple, Ivan der Tak et all (|2002j ) detected triply-deuterated 
NH3 toward IRAS 4A in a 25" beam, where we find 
Rd = 0.16. 

Ten of the sources analysed by lEmprechtinger et al.l 
overlap with protostellar cores in our Perseus sample. 
A direct comparison of the deuterium fractionation val- 
ues for the overlapping targets shows that their results 
are greater, on average, than those reported here, with 
varyin g discrepancies that range from a factor of 0.8 
to 3. lEmprechtinger et al.l derive both A^(N2lI+) and 
A^(N2D+) from observations taken with the IRAM 30 m 
telescope based on the N2II+ J = 1 — and N2D+ 
J = 2 — 1 transitions, respectively. While the N2II+ 
observations (and consequently the derived column den- 
sities) thus match ours in transition and beam size, the 
angular resolution of the N2D+ 2-1 observations (16") is 
significantly better than both our N2D+ 2-1 (41") and 
3-2 (30") observation s. In all but one overlapping source, 
lEmprechtinger et al.l report greater integrated N2D+ 2- 
1 line intensities than found in our observations. This 
suggests that the filling factor of the N2D+ emission, 
here assumed to be unity, is < 1 for some sources. We 
discuss the filling factor assumption and its effect on our 
analysis further in ij4.1.2l 



3.2.1. Trends in the deutenum fractionation 

We plot in Figure [7] the N2II+ deuterium fraction 
against several core parameters thought to be indicators 
of core evolutionary stage: the N2H+ 1-0 velocity dis- 
persion, ffy, the kinetic gas temperature, Tk, and the 
column density of H2, Nj Ki). Where targets coincide 
with a SCUBA core in the iKirk et al.l (|2006) catalogue 
(51/64 sources), we also plot Rd against the core central 
density, Uc, from a Bonnor-Ebert (BE) fit to the data. 
For all parameters plotted, we test for significan t correla- 
tions using Spearman's p rank correlation test (jConoveii 
[T999I) . 

The H2 column density, A^(Il2), is calculated from 
SCUBA 850 liui continuum emission (S^; described in 
IKirk et aI]|2006D . following the standard modified black- 
body equation 7V(H2) = I {ftmlJ''mHK^B^Tci). Here, 
fi = 2.33 is the mean molecular weight, mn is the mass 
of hydrogen, = 0.02cm^g~^ is the dust opacity per 
unit mass at 850 /im, and B^{Td) is the Planck function 
at the dust temperature Td. We sum flux over a 30" 
beam {Qm) for each targ et. The dust opacity i s taken 
from the dust models of lOssenkopf fc HenningI (|1994D 
for agglomerated grains with ice mantles (0II5). This 
value is consistent with observational results in low mass 
star forming cores (Shir lev et al.l l2002). We set the dust 
temperature equal to th e NH.'^-derived gas temp erature 
where cores overlap with iRosolowskv et al.l ()2008D . Else- 
where, we set Td = 11 K, the mean of the Perseus sample. 
Our iV(H2) values thus differ slightly from those reported 
for overlapping cores by iJohnstone et al.l (|2010( ) , who set 
Td — 11 K (the mean value for all cores in their study) 
for all cores. There are several sources of uncertainty in 
the A^(Il2) values. First, the calibration uncertainty in 
the continuum fiux values is ^ 20%. Second, the dust 
properties are not known with high precisi on. For ex- 
ample, estimates of ksso can vary by ~ 3 ([Shirlev et al.l 
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Fig. 5. — Top: Distribution of Rn values found for Perseus cores. 
The total distribution is shown in black; blue and red histograms 
show the distribution for starless and protostellar cores, respec- 
tively. Mean values are shown with arrows. Bottom: Total Ro 
distribution as in (top) , with blue and red histograms showing the 
distribution for isolated and clustered cores, respectively. Mean 
values are shown with arrows. 



[2000h . while Td may differ from Tk- At the densities typ- 
ical of the cores studied here (n > 10^ cm^*^ ), however, 
we ex pect that the dust and gas are coupled ()Goldsmithl 
1200 It) . The BE analysis independently derived core tem- 
peratures of T = 10 — 19 K for the Perseus cores, in gen- 
eral agreement with the derived Tk values we use here. 
In addition, while ksso may not be known to high preci- 
sion, we do not expect large variations in Kgso between 
the cores given their similar temperatures and densities. 
We estimate the H2 column densities have uncertainties 
of factors of a few. 

We find no trend in the deuterium fractionation of 
the cores as a function of N2H"'" line width ( Figur e 
[TK- Sim ilar results were fo und b y iCrapsi et all (pOOSh . 
and by [Roberts fc Millaij ()2007[ ) toward a sample of 
low- and high-mass protostellar cores. Towar d cores 
in massive star forming regions, iFontani et al.l (1201 If ) 
find a faint anticorrelation between Rn and N^Hj line 
width (J = 3 — 2 transition). iChen et al.l (|2011t ) also 
find an anticorrelation between Ru and the N2D"'^ line 
width (J = 3 — 2 transition) toward massive proto- 
stellar/cluster cores. Anticorrelation between deuterium 
fractionation and line width may be expected if the 
turbulence typical of the interstellar medium has dis- 



sipated within the most evolved dense cores, resulting 
in gas motions that are nearly purely thermal. Dense 
gas within gravitationally unstable cores may also, how- 
ever, have significant infall motions (^0.1 km s~^) trace d 
byN2H+ (iSchnee et all [20071 IPTFrancesco et atl lMT). 
which would add a broader component to the observed 
line widths and may confuse any trend with deuterium 
fraction. 

We additionally find no trend in the deuterium frac- 
tionation of the cores as a function of gas tempera- 
ture, Tk (Figure Uh)- Multiple studies have found 
significant anticorrelation between temperatu r e an d 
deuterium fractionat i on, i n cluding iChen et al.l (|20Tll) ; 
lEmprechtinger et al.l (|2009D : IFontani et al.l (|20ill) . This 
anticorrelation is expected for several reasons. First, the 
reaction that forms H2D+ from is exothermic, and 
at low temperatures can only proceed in the forward di- 
rection. This results in a reservoir of H2D+ to react 
with molecules in the gas phase, and produce deuterated 
species. At temperatures greater than ~ 25 K, H2D+ 
is able to react with H2 and th e H2D+ abundance de- 
creases ([Roberts fc Millaill2007[ ). Second, at low temper- 
atures (~ 20 — 30 K) CO becomes depleted from the gas 
phase through freeze-out onto dust grains, removing a de- 
struction pathway for H2D+ (and N2D+) and thus allow- 
ing enhanced abundances of deuterated species. When 
a protostar forms and begins warming the nearby gas, 
CO can desorb from the dust grain, returning to the 
gas phase and again reacting with H2D+ (and N2D+) to 
suppress deuterium fractionation in molecular species. In 
the studies already mentioned, the decrease in deuterium 
fractionation is largely seen at temperatures > 16 — 25 K. 
Because all the Perseus cores studied here have gas tem- 
peratures Tk < 17 K, it is reasonable that Figure [7b 
shows no strong correlation between deuterium fraction- 
ation and temperature. 

We note, however, that the measured gas temper- 
atures for our targ e ts wh ich overlap with those of 
lEmprechtinger et al.l (|2009[ ) are significantly colder (of- 
ten by factors of 2 or more) than the dust temperatures 
the authors determined through spectral energy distribu- 
tion (SED) fitting, and where the authors find a strong 
anti-correlation between Ro and T^. At the densities 
typical of the cores studied here {n > 10^_cm~fV we ex- 
pect that the dust and gas are coupled (|Goldsmith| [2001). 
Since warm dust radiates more efficiently, however, it is 
possible that the average derived dust temperature of a 
core with a warm central component (i.e., due to heat- 
ing from a young protostar) will be greater than that 
probed by NH3 observations. Compared with our results 
at larger angular resolution, this suggests that the aver- 
age deuterium fraction of dense gas-tracing molecules in 
the protostellar cores is dominated by the physical con- 
ditions close to the protostar. A s imilar result was found 
previously bv lFriesen et al.l ([2O1O0 toward the Ophiuchus 
B core, where the deuterium fraction in N2II+ and II2D+ 
abundance both decrease with increasing proximity to 
protostars embedded within the core. 

The large discrepancy between Td and Tk could also 
be due to the difficulty in fitting dust SEDs with few data 
points, however, where anticorrelations between Td and 
the dust opacity power la w index, /3, ca,n be rn imicked by 
noise m the d ata po ints ([Shettv et al.ll2009al[bD . While 
lEmprechtinger et al.l model the SEDs of their sources 



12 



Friesen et al. 




48m00s 3h47m30s 44m00s 3h43m30s 




28mOOs 3h27m30s 3h25m30s 

RA (J2000) 



Fig. 6. — 850 iim continuum emission toward well-known star-forming regions within the Perseus molecular cloud taken with SCUBA 
on the JCMT ijKirk et al.ll2009 l. In Perseus, the regions shown run approximately east to west from top left to bottom right. To better 
show the cloud structure, the greyscale has been truncated at 0.2 Jy beam~^, while peak values reach ~ 8 Jy beam~^ toward several regions 
in NGC 1333. Grey contours show flux densities of 0.1 Jybeam"^ and 0.5Jybeam~^. Overlaid are circles (30" diameter) showing the 
positions of 54 of 64 cores observed in this study, with colours indicating their derived N2H+ deuterium fraction, Rr>. Cores with Ro > 0.1 
are shown in red, i?£) < 0.1 are shown in orange, and N2D+ non-detections are shown in blue. The remaining ten cores are located outside 
these well-known regions: one shows Rjj > 0.1, six show i?j5 < 0.1, and three do not have N2D+ detections. 
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using /3 = 2, consistent with the 0H5 dust opacity model 
described above including a strong gradient in Td in their 
models, they then determine a single dust temperature 
Tci from an isothermal fit to the data, assuming /3 = 1. 
The higher dust temperatures determined through this 
analysis may better reflect the physical conditions of the 
dust close to the protostar than the Tk measurements. 
Since our observations provide a beam-averaged measure 
of the deuterium fraction, however, we conclude that the 
Tk measurements provide the best estimate of the core 
temperature for our analysis. 

We find in Figure [7h that the N2H"'" deuterium frac- 
tion tends to increase with the column density of H2, 
7V(H2). Although there is substantial scatter, a statis- 
tical test shows the correlation is significant (Spearman 
parameter p = 0.29, with a significance of 0.04, where 
low values indicate greater significance) . Considering the 
starless and protostellar populations separately, we find 
that the deuterium fractionation in the protostellar pop- 
ulation is correlated with A^(H2) (p — 0.36, significance 
0.04 omitting the embedded sources, and p — 0.27, sig- 
nificance 0.12 including the embedded sources). There 
is no significant correlation in the starless population. 
Similarly, a general increase in Ro with central density 
Uc is seen in Figure [7Ji when considering the entire core 
sample, where p = 0.23 with a significance of 0.1. In this 
case, however, the correlation is stronger when analysing 
the starless cores separately (p = 0.55, significance of 
0.02). 

3.2.2. Deuterium fractionation and CO depletion 

Models of dense gas chemistry show that the de- 
pletion of CO from the gas phase increases the rate 
of deuterium fractionat ion reactions in the gas phase 
(jPalgarno fc Lepij|1984[ ). At low temperatures, the deu- 
terium fractionation rate is stron gly de pendent on the 
CO depletion (|Casem et al.l [200ll2008h . Recent stud- 
ies have indeed found strong correlations between the 
deuterium fractionation of N2H"'" (and other molecu- 
lar species) and the amount of CO depletio n, /d, in a 



star-f orming core (fCa scUi ct al. 2002; Crapsi et al.l 



2004 . 



2005; Casclh et al.l [2008; Emprcc htinger et al.l 12009 : 
Chen et al. 201ll) . Here wo determine the beam-averaged 
CO depletion tow ard the Perseus c ores using the C^^O 
data presented in IKirk et all (j2007D . The CO depletion 
is then fo = x,a„(Ci80)/a;(Ci80), where XcaniC^^'O) 
is the canonical abundance of C^^O relat ive to H2 in 
the IS M. To best compare our res ults with iCrapsi et al.l 
(|2005f) . we set x^an = 1-7 x 10""^ (|Frerking et al.lll982D . 
but note t hat other values of Xcan exist in the literature 
(e.g.. ILacv ct al. 1994). We discuss further the effects of 
using the Frcrking et ahl value in fS 



To calculate iV(C^^O), we assume that the C^^O emis- 
sion is optically thin, and set Tex = Tk for all sources. 
We then derive 7V(C"'^*0) in the optically thin case us- 
ing th e line integrated intensity, following iCaselli et al.l 
(|2OO20, and determine Jd for each source using the pre- 
viously derived iV(H2) ( §3.2.ip . We expect the assump- 
tion that Tex = Tk to be reasonable for the Perseus cores, 
given that C^^O traces lower density gas than NII3. At 
lower densities, the gas-dust coupling is less efficient, and 
the gas kinetic temperature may be elevated above the 
dust temperature by photoelectric heating. In this case, 



the C-'^^O kinetic temperature, Tk,cOi may be somewhat 
higher than for NH3. Since the C^^O is likely subther- 
mally populated, however, T^x < Tk.co- If the C^^O 
emission has moderate opacity, iV(C^^O) would be un- 
derestimated (and /d would be overestimated) by ~ 30% 
for r = 0.5, and ^ 60% for r = 1. The uncertainty in the 
CO depletion also depends linearly on the uncertainty in 
iV(Il2), which may be relatively large as described above. 
For example, if we have overestimated the dust emissiv- 
ity at 850 pm for the Perseus molecular cloud, our de- 
rived 7V(H2) values will be systematically low, as will 
our derived CO depletion factors. Variations in the dust 
properties between the cores will introduce scatter in the 
results. Two cores lie off the area mapped by SCUBA, 
allowing no measurement of A^(H2), and consequently 
Id. 

On average, Perseus cores are characterized by low 
values of fo- CO depletion factors fa > 5 are found 
only toward three protostellar cores, including NGC 1333 
IRAS 4A and IRAS 4B (our sources 106 and 103), where 
/n = 8 and fo = 19, respectively. Both IRAS 4A and 
IRAS 4B h ave been previousl y shown to have high deple- 
tion of CO (IBlake et al.llT995h . in agreement with our re- 
sults. More recently, by mod eling the observed emission 
from seven C^^O transitions. lYildiz et all (|2012D showed 
that CO abundance profiles toward both sources must 
include a factor of > 10 drop in abundance between 
~ 5 X lO^AU and - 10"* AU (where lO^AU is ~ 40" 
at the distance of Perseus). As expected, IRAS 4A and 
4B are also two of the Perseus cores that show the high- 
est N2II"'" deuterium fraction levels, with Ru = 0.16 and 
Rd = 0.20, respectively. 

All other cores have CO depletion factors /d ^ 5, with 
a mean fu = 2.5 for the sample. We additionally find 
/d < 1 for 12 cores, meaning that these cores contain 
more C^^O relative to H2 than expected from the canon- 
ical value. Seven of these over-abundant C^®0 cores lie 
in the eastern part of Perseus. For example, cores in the 
B5 and IC348 regions have a mean = 1.0 and 1.3, re- 
spectively, with small standard deviations of 0.5 and 0.7. 
Cores in NGC 1333, L1448, L1455 and Bl show greater 
average fo values, with a mean /£> 3 — 4, and sub- 
stantially larger variations in fjj between cores in each 
region. The results are suggestive of a gradient in the 
C^®0 depletion across Perseus, but we note that cores in 
B5 and IC348 together include only 13 of our sample of 
64 cores. 

We plot in Figure [8] the N2H+ deuterium fraction 
against the CO depletion factor for each core. In gen- 
eral, the Figure shows that Ru tends to increase with fu, 
although the trend is most significant for the protostel- 
lar sample, omitting the most embedded sources (Spear- 
man's p = 0.34, with a significance of 0.05). On average, 
protostellar cores are also characterised by significantly 
greater CO depletion (mean — 3.2) than starless cores 
(mean fjj = 1.2), with a < 1% KS2 probability that 
the distributions are dr awn from the same sample . For 
ten cores in our sample. FEmprechtinger et al.l ()2009l ) also 
calculated the CO depletion factor using observations of 
C^^O 1-0 (22" HPBW). They find fo values that are 
generally within factors of ~ 2 of the results presented 
here. Most (but no t all) of our fu values are greater than 
those presented by lEmprechtinger et al.l This is reason- 
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Fig. 7. — Trends in the N2H+ deuterium fractionation, Rd, with N2H"'" hne width (uv), gas temperature (Tk), H2 column density 
(Af(H2)), and the log of the central density, ric- The gas temperature is calculated from analysis of NH3 (1,1) and (2,2) inversion transition 
emission fRosolowskv et al. 2008). Af(H2) is calculated from 850 /xm continuum measurements IIKirk et al.ll2009 ). recalculated here using 
the NHs-dorivod gas temperature. The central density is derived from a Bonnor-Ebert sphere analysis of the 850 fim continuum emission 
toward the cores (iKirk et al.iL200ti,) . In all plots, starless and protostellar cores are shown by black and grey (colored red in the online 
version) triangles, respectively, while deeply embedded sources are shown as open triangles (colored blue in the online version). 



able considering our C^^O observations have finer angu- 
lar resolution, and are thus sampling emission closer to 
the core center, encompassing less of the lower density 
gas where we expect lower CO depletion. 

It is interesting that many of our cores show low fu 
values. Low CO depletion will stron gly decrease the deu - 
terium fractionation, and models bv lCaselli et al.l ()2008f ) 
predict Rd < 0.3 for fu < 10, in agreement with our 
results. Trends between Ru and fo found in previous 
studies ofte n show a marke d increase in Rp at fo > 10 
(jCrapsi et al. 2005; Emprc chtinger et al.|[2009[ ). At low 
fn < 10. [Chen et al.i (|2011l ) find a weak increasing trend 
in Rjj with /u in their study of mas sive protostellar 
cores. Recently, iFord &: ShirlevI (j2011f) show, through 
radiative transfer modeling of C^^O data toward eight 
starless cores in Taurus, that fjj is often very poorly 



constrained, largely due to uncertainties in the canonical 
abundance of C^^O and in the core temperature. They 
find that it is difficult to distinguish between — 10 
and fo — 1000 in their models. In Perseus, however, 
even if the absolute fo values are uncertain, the finding 
that Perseus cores have low CO depletion on average is 
likely robust, with the caveat that uncertainties in the 
dust emissivity may have a systematic effect on our re- 
sults. We discuss the implications of the low depletion 
values in Perseus in more detail in m.2\ 

3.2.3. Deuterium fractionation and protostellar properties 

Finally, we examine the N2II+ deuterium fraction of 
the protostellar Perseus cores with respect to the prop- 
erties of the embedded young stellar ob jects (YSOs). In 
their analysis of YSOs in Perseus, lEnoch et al.l (|2009( ) 
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Fig. 8. — The N2H+ deuterium fractionation for Perseus cores as 
a function of CO depletion, fo, as measured tlirougli observations 
of C^*0 2-1 emission at 11" resolution. Starless and protostel- 
lar cores are shown by black and grey (colored red in the online 
version) triangles, respectively, while deeply embedded sources are 
shown as open triangles (colored blue in the online version). 



determine the bolometric temperature (Ti,oi) and lumi- 
nosity {Li,oi), as well as envelope mass (Mem,), based 
on spectral energy distribution fits to 2MASS, Spitzer 
IRAC and MIPS photometry, and Bolocam 1.1mm data 
{Fi.imm)- Younger sources are expected to have lower 
luminosities and temperatures, greater remaining enve- 
lope mass, and low ratios of Lboi to submillimetre flux 
density. Of the 43 protostellar cores surveyed here, 39 
have associated Tboi, Lboi, M^nv and i^i.imm measure- 
ments, of which we detect 32 in N2D+ emission. Within 
this sample, 21 are identified as Class sources and 18 
as Class I, whe re Class sou rces are defined as sources 
with Tboi < 70 (jEnoch et al.t ). We show in Figure |9] the 
deuterium fraction of these 39 sources as a function of 
Lboh Tboh Menv and Lboi/Fi,imm- Note that the three 
deeply embedded protostellar sources, which were not 
identified in the infrared surveys, have estimated internal 
luminosity u pper limits Ljnt S, - 01 Lp) (Per-Bolo 45 and 
Per-Bolo 58; iSchnee et all [20121: lEnoch et ahllMol ). or 
bolome tric lum inosity upper limit Lboi < 0.1 £© (L1448- 
IRS2E; [Chenet al. 2010). 

While the mean Rd for the Class sources is slightly 
greater than that for the Class I sources (Rd = 0.08 and 
0.07, respectively), the difference is not significant. We 
find no significant trends in Rd with protostellar Lboi-, 
Menv, or Lboi/ Fi,imrm but do find an anti correlation be- 
tween Rd and Tboi (Spearman p = —0.25, significance of 
0.07). In addition, the sources with greater Rd values are 
found toward sources with Menv ^ 0.8 M©, despite the 
lack of a clear correlation between Rd and Menv ■ Fur- 
thermore, all the protostellar sources in the lEnoch et al.l 
sample that we did not detect in N2D"'" emission have 
higher Tboi values, identifying them as older. Class I 
sources, while six of seven undetected sources have 
Menv £ 0.5 M0 (see Rd upper limits in Figure [nt). Fi- 
nally, sources with high Lboi/ Fi,imm values tend to have 
the lowest Rd measurements, although there is signif- 
icant scatter at Lboi/ Fi imm ^ 2 — 3L0 Jy~^. Taken 



together, these plots show that the deuterium fraction 
in N2H+ in the protostellar Perseus cores is consistent 
with the expectation that Rd will drop as the protostar 
warms its immediate environment (releasing CO from 
dust grains into the gas) and additionally loses envelope 
mass to accret ion and outfiow, and i s in ag reement with 
the findings of lEmprechtinger et al.l ()2009l ). 

4. DISCUSSION 

4.1. Potential complications of a pointed survey 

4.1.1. Excitation effects 

In §3.2) we set the excitation temperature of the N2D+ 
emission equal to that of the N2H+ emission for each 
core to determine iV(N2D+). This assumption has been 
shown to be accurate in several dense cores through high 
sensi tivity observations (jCaselli et al.l l2002t iGerin et al.l 
[2OOI . For those sources in this study where we detect 
both N2D''" transitions, we can instead simultaneously 
determine both the column density of N2D+, A^(N2D+), 
and Tex , under the assumptions that both lines are opti- 
cally thin and emission fills the beams, and requiring the 
excitation temperatures and total N2D"*' column density 
derived from each transition to be equal. Uncertainties 
in the derived properties were calculated using the for- 
mal uncertainties in the measured integrated intensities 
in a Monte Carlo analysis. In this analysis, using the 
column density calc ulation for optically thin emission in 
ICaselli et al.l (j20 02f). the excitation temperature derived 
from the two species is solely dependent on the ratio of 
the integrated intensities, W2/W1, where 



W2 _ 9u.2Xl A2 J,,{Tex) 

Wi gu.iXl Al J,,{Tex) 



J U2 (Fbi 



Ju2 {Fbg ) 



1 - exp{-hi'2/kTex) exp{-Ei^2/kTex) 
1 - exp{-hiyi/kTex) exp{-Eii/kTex)' 



(3) 



Here, f/„_i and gu,2 are the upper state statistical weights 
for the two transitions, Ai and A2 are the transition wave- 
lengths, Al and A2 are the Einstein A coefficients, Ei^i 
and Ei^2 are the lower state energy levels. Ji,-^{Tex) and 
(Tex) are the equivalent Raleigh- Jeans excitation tem- 
peratures for the two transitions at excitation tempera- 
ture Tex', J VI (Tbg) and (Tbg) are the same for the back- 
ground temperature, Tbg = 2.73 K. The derived N2D+ 
column density depends both on Tex and the measured 
integrated intensities. 

The resulting N2D+ Tex values are often, but not 
always, in agreement within uncertai nties with those 
measured from N2H"'' 1-0 observations (jJohnstone et al.l 
l2010f l. Further separation of these sources into protostel- 
lar and starless populations reveals that the N2D+ Tex 
values for protostellar sources are greater, with a mean 
Tex = 8.4 K and a standard deviation of 2.5 K, com- 
pared with starless sources, where the mean Tex = 6.7K 
with a standard deviation of 1.8 K. There is little dif- 
ference, however, in the N2II+ Tex between the starless 
(mean Tex — 5.9 K) and protostellar (mean Tex — 6.0 K) 
populations, with the N2H+ and N2D+ excitation tem- 
peratures consistent within uncertainties for the starless 
cores. Greater N2D+ excitation temperatures for proto- 
stellar cores, if real, suggest that the N2D+ observations 
may trace somewhat warmer or higher density gas near 



16 



Friesen et al. 




Fig. 9. — The N2H+ deuterium fractionation of protostellar cores in Perseus against (a) Lfcof, (b) Ttoi, (c) the envelope mass Menv, and 
(d) the rat io of Lf,ai to t otal 1.1mm flux. L(,o;, T^oj and M^nv values are taken from lEnoch et ahl (2009), wrhilc the total 1.1mm flux is 
taken from [Enoch et al] 1200 6). These catalogs do not include recently discovered very embedded sources. In all plots, Class and Class 
I sources, as determined by Enoch ot al., are shown grey (colored red in the online version) and black points, respectively. For clarity, we 
show one-sided x-axis error bars on the Ljoj plots. 



the embedded protostars than N2H+ 1-0. Ahernatively, 
however, this increase in measured T^x may be caused in- 
stead by variation in the emission fiUing factor within the 
telescope beams for N2D+ 2-1 and N2D+ 3-2 emission. 

4.1.2. Source sizes 

The excitation analysis is strongly dependent on the 
assumption that the emission for both species fills the 
respective beams (the filling fraction for both transitions 
is unity). The FWHM beam widths are 41" and 30" for 
the N2D"'" 2-1 and N2D+ 3-2 transitions, respectively. 
To test this, we included a filling factor / equal to the 
ratio of the beam areas between the N2D+ 3-2 and N2D+ 
2-1 observations and determined again T^^x and A^(N2D+) 
following Equation |31 Essentially, this calculation as- 
sumes that the size of the emitting region is similar to 
the N2D+ 3-2 beam FWHM. When incorporating the 



beam ratio in the analysis, we find the excitation tem- 
perature of the N2D+ emission is within uncertainties 
of the N2H+ Tex for most sources. This agreement sug- 
gests that the size of the N2D+ emitting region is indeed 
smaller than the N2D+ 2-1 beam, and likely similar to 
the N2D+ 3-2 beam. 

Furthermore, if we assume the excitation of the N2D+ 
emission is equal to the N2H+ excitation, we can de- 
termine the relative filling fraction of the beam between 
the two N2D+ transitions, although we cannot deter- 
mine accurately the size of the emitting region if it 
is significantly smaller than 30". The resulting N2D+ 
source sizes range from ~ 25" to ~ 40", with a mean 
of 28". This is in broad agreement with the source sizes 
from the SCUBA 850 /im continuum measurements, al- 
though individual source sizes do not agree well. Previ- 
ous studies have shown that the emitting region of N2D+ 
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can be smaller tha n the SCUBA source size (see, e.g., 
iCaselli et al.l 120021 ). This result is expected given that 
N2D+ should be emitted from regions of high density, 
whereas the dust continuum emission will have some con- 
tribution from lower density material (although submil- 
limeter and millimeter continuum observation methods 
generally result in the loss of emission on large angular 
scales, removing some of the lower density component). 
In addition, we are able to directly compare the observed 
integrated intensity of N2D"'" 2-1 toward seven proto- 
stellar sources in our samp le which overlap with IRAM 
N2D+ 2-1 observations bv lEmprechtingeT et al.l pOOO). 
The measured N2D+ integrated intensities agreed for 
only one source (Per4N in Emprech tinger et al.; target 97 
in our catalogue). For the rest, the lEmprechtinger et al.l 
()2009f ) integrated intensities were significantly greater. 
The resulting filling fractions for the overlapping targets 
in N2D'*' 2-1 emission range between 0.23 and unity, or 
source sizes of 20" to at least 41" (the FWHM of our 
N2D+ 2-1 observations). 

We find that the emitting region sizes derived from 
these tests are significantly larger for protostellar cores 
than for starless cores, and we additionally find a pos- 
itive trend between source size and core concentra- 
tion. This is in agreement with previous results that 
the mo re massive cores i n Per seus tend to be proto- 
stellar (jHatchell fc Fulled 120081 ). and that highly con- 
centrated cores are also more likely to be protostellar 
(jWalawender et al.l [2005^. In sum, protostellar cores are 
more massive, more concentrated, and show N2D+ emis- 
sion over larger areas, but on average the deuterium frac- 
tionation in N2H-I- remains the same as in starless cores. 

For most of the cores, however, we derive Rd from 
the integrated N2D+ 3-2 intensity. The effect of vari- 
able filling fraction is therefore lessened, where the small- 
est filling fraction of 0.23 calculated for N2D+ 2-1 from 
the lEmprechtinger et al.l results increases to 0.44 in the 
N2D+ 3-2 beam. The agreement between the N2H+ and 
N2D"'" Tex values when incorporating the beam ratios in 
Equation [3] suggests that in most cases, the N2D+ 3-2 
beam is similar in size to the emitting regions. We ex- 
pect, therefore, that a low filling factor would decrease 
our derived Rd values by a maximum of ~ 2 for a few 
sources. While we therefore expect that the introduced 
scatter should be small, it may be significant if the N2D+ 
source size is correlated with other core properties. This 
effect could be probed through high resolution maps of 
the protostellar sources in N2II+ and N2D+ emission. 

4.1.3. Offsets of N2D^ emission from the dust continuum 
peak and abundance gradients 

Offsets between the position of peak N2D+ emission 
and our pointing locations would be an additional source 
of scatter in our results. As discussed in H2.ll our target 
list is dominated by cores identified based on SCUBA 
850 /xm continuum emission. Through additional obser- 
vations in a cross pattern around many of their targets, 
iKirk et al.l (|2007| ) showed that the peak N2H+ emission 
toward these cores is coincident with the continuum peak 
for nearly all their SCUBA sources. The targets se- 
lected based on extinction data were more likely to have 
brighter N2H+ emission offset from their central point- 
ing location. Of the 13 sources in our sample not se- 
lected from SCUBA data, only three show substantially 



stronger N2H+ emission ^ 25" offset to our pointing lo- 
cations. We thus are confident we are probing the peak of 
the N2H"'" emission for the large majority of our sources. 

Some studies have shown that the emission of deuter- 
ated species can be significantly offset from associ- 
ated submillimeter continuum emis sion, or even frorn . 
their non-deuterat ed counterparts (jPagani et al.l 120051 : 
iFriesen et al.ll2010l ). Given the beam size of the obser- 
vations presented here (30" and 41"), the N2D+ emis- 
sion would have to be substantially offset from the N2II+ 
emission to be problematic for this analysis. Based on 
850 /im continuum maps, the median effective radius of 
our targets is 28". Since N2D+ emission should remain 
associated with the core, even if offset from the central 
peak as measured by submillimeter continuum or N2II"'" 
emission, we expect scatter due to emission offsets to be 
small. Small maps of N2 11+ and N2D+ toward the cores 
would be necessary to quantify this effect. 

4.2. Comparison with chemistry models 

Despite the lack of clear trends in the deuterium frac- 
tionation of the Perseus cores with many parameters 
that may trace an evolutionary sequence, the data nev- 
ertheless inhabit a parameter space that needs to be ac- 
counted for in models that combine physical and chem- 
ical evolution of low-mass star forming cores. In gen- 
eral, our results are in broad agreement with published 
steady-state and dynamic chemical models of core evolu- 
tion. The deuterium fractionation values for the Perseus 
co res are bracketed by the slow and fast collapse models 
of lAikawa et al.l (|2005l ). who show that cores that col- 
lapse quickly will show less deuterium fractionation in 
N2H+ compared with cores that collapse more slowly. 
The radial velocity profiles predicted for the fast col- 
lapse models are large, however, compared with typi- 
cal observational results. While large infall speeds have 
been found for the protostellar sources NGC 1333 IRAS 
4A and 4B (0.7kms~^ and 0.5kms~^, respectively; 
iDi Francesco et all 120011 ). mo re typic al infall speeds are 
~ 0.05kms-i to O.lkms-i (|Ma.rdon cs 1998). Steady- 
state chemistry core models can reproduce the range of 
Rd values fo und in Perseus give n variations in the model 
parameters (jCaselli et al.|[2008() . The dominant param- 
eters that affect the deuterium fractionation at the low 
temperatures (T < 17 K) characteristic of our targets 
in Perseus are the CO depletion, the core density, the 
fractional abundance of dust grains with respect to H2, 
the cosmic-ray ionizati on rate, and the ratio of ortho-Hg 
to par a-H9 in the gas (jPagani et al.lll992HGerlich et al.l 
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'Casel li et al. 1 120081: IRoueff et al. 1 120051 IRoberts fc Millari 
2000). 

The amount of CO depletion in cold, dense gas strongly 
impacts the level of deuterium fractionation that ca n 
be reached (jDalgarno fc Lepp|ll984l: iCaselh et al .1120081) . 
Using a simple chemical model. iKeto &: Casellil ([2008f ) 
calculate the expected steady-state CO abundance as 
a function of density and temperature in starless cores. 
Given their average densities and temperatures, we ex- 
pect CO depiction factors /d '--^ 6 — 8 for the Perseus 
cores. In ij3.2.2[ however, we find that the CO deple- 
tion factors for most of the Perseus cores are low (mean 
/d = 3.8), with all but three cores having depletion fac- 
tors Jd < 5. In fact, we find an over-abundance of CO 
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TABLE 3 

Deuterium fractionation and CO depletion in the Perseus cores 



oource 


iV (^iN2tl ) 


i V (IN 2 ; 




iV(H2) 


iV \ \J ) 


ID 




X 10 cm 


XlU cm 




Xl02^cm-2 


1 nl4 — 2 

X 10 cm 








Dstsctsd in 


both N2D+ 2-1 and 3-2 








97 


19.2 (1.4) 


0.198 (0.042) 


1.7 


20.9 


1.41 


21 


81 


6.2 (0.3) 


0.077 (0.016) 


1.6 


29.9 


0.92 


23 


56 


7.8 (1.5) 


0.140 (0.038) 


1.3 


36.1 


0.63 


25 


125 


23.8 (0.7) 


0.190 (0.038) 


4.6 


37.3 


2.10 


26 


105 


8.2 (0.3) 


0.078 (0.016) 


2.4 


17.1 


2.38 


27 


42 


4.9 (0.8) 


0.116 (0.029) 


1.1 


25.9 


0.73 


36 


59 


2.6 (1.0) 


0.044 (0.019) 


0.3 


8.4 


0.66 


71 


355 


33.9 (1.1) 


0.096 (0.019) 


6.7 


24.6 


4.65 


72 


248 


17.1 (0.5) 


0.069 (0.014) 


6.5 


41.6 


2.64 


73 


343 


31.7 (1.2) 


0.092 (0.019) 


4.0 


32.7 


2.09 


74 


128 


12.1 (1.3) 


0.095 (0.022) 


2.0 


14.6 


2.28 


75 


203 


9.0 (1.9) 


0.044 (0.013) 


1.8 


22.2 


1.35 


77 


115 


7.4 (1.5) 


0.064 (0.018) 


1.3 


18.2 


1.23 


86 


171 


14.2 (0.6) 


0.083 (0.017) 


6.3 


22.7 


4.72 


90 


66 


10.7 (0.5) 


0.163 (0.033) 


3.0 


13.7 


3.74 


92 


193 


11.3 (1.5) 


0.059 (0.014) 


0.8 


7.4 


1.89 


95 


126 


9.5 (0.8) 


0.075 (0.016) 


1.7 


14.9 


1.98 


96 


135 


13.8 (1.2) 


0.102 (0.022) 


1.3 


11.4 


1.96 


97 


116 


7.0 (0.6) 


0.060 (0.013) 


1.7 


12.7 


2.23 


99 


48 


5.3 (0.7) 


0.111 (0.026) 


1.3 


11.3 


1.98 


100 


80 


5.2 (0.5) 


0.065 (0.014) 


1.3 


31.1 


0.71 


102 


229 


14.6 (0.7) 


0.064 (0.013) 


2.2 


35.2 


1.06 


103 


104 


20 3 (n 6") 


0.196 (0.039) 


12.5 


11.3 


18.89 


104 


142 


4.6 (0.3) 


0.032 (0.007) 


4.1 


38.8 


1.81 


106 


147 


22 Q (n 8") 


0.156 (0.032) 


18.4 


38.6 


8.08 


107 


209 


24.4 (0.5) 


0.117 (0.024) 


3.6 


27.0 


2.28 


109 


77 


9.2 (0.5) 


0.120 (0.025) 


1.5 


36.3 


0.72 


110 


150 


10.8 (0.3) 


0.072 (0.015) 


3.0 


45.6 


1.11 


111 


233 


6.3 (0.3) 


0.027 (0.006) 


1.6 


16.5 


1.68 


112 


233 


7.7 (0.3) 


0.033 (0.007) 


7.9 


56.4 


2.39 


113 


176 


14.3 (0.5) 


0.081 (0.016) 


3.3 


74.9 


0.75 


115 


128 


5.4 (0.5) 


0.042 (0.009) 


0.6 


19.0 


0.57 


116 


113 


7.6 (0.3) 


0.067 (0.014) 


2.8 


48.7 


0.98 


118 


161 


6.0 (0.5) 


0.037 (0.008) 


4.5 


39.7 


1.91 


121 


100 


5.2 (0.6) 


0.052 (0.012) 


1.3 


15.1 


1.49 


122 


210 


7.3 (0.7) 


0.035 (0.008) 


2.3 


14.7 


2.68 


123 


208 


8.8 (1.0) 


0.042 (0.010) 


2.3 


21.4 


1.80 


124 


124 


6.3 (0.8) 


0.051 (0.012) 


1.4 


24.0 


1.01 


134 


158 


7.0 (0.6) 


0.044 (0.010) 


2.3 


10.1 


3.81 


135 


206 


9.0 (1.1) 


0.044 (0.010) 


2.0 


6.5 


5.26 


136 


89 


3.3 (1.1) 


0.038 (0.014) 


1.1 


7.3 


2.57 


146 


132 


13.4 (1.1) 


0.101 (0.022) 


3.0 


22.1 


2.29 


148 


169 


22.8 (1.1) 


0.135 (0.028) 


4.0 


9.5 


7.07 


149 


332 


25.1 (0.7) 


0.076 (0.015) 


10.5 


42.6 


4.20 


150 


253 


15.8 (1.8) 


0.063 (0.014) 


2.5 


17.6 


2.41 


152 


313 


10.4 (1.5) 


0.033 (0.008) 


4.3 


24.5 


2.95 


Detected in N2D+ 2-1 only 


36 


59 


1.7 (0.3) 


0.029 (0.008) 


0.3 


8.4 


0.66 


147 


83 


5.0 (0.4) 


0.060 (0.013) 


0.5 


22.3 


0.37 


Detected in N2D+ 3-2 only 


3 


145 


4.1 (0.7) 


0.029 (0.008) 


0.3 


13.1 


0.43 


4 


145 


4.1 (0.5) 


0.028 (0.007) 


1.6 


23.3 


1.19 


15 


55 


4.4 (0.9) 


0.080 (0.023) 


1.2 


23.3 


0.87 


19 


66 


1.6 (0.4) 


0.024 (0.008) 


1.5 


19.1 


1.35 


143 


114 


5.0 (1.4) 


0.044 (0.015) 


1.0 


5.3 


3.25 






iNon-dctcction upper 


limits 






5 


33 


< 1.5 


< 0.030 


0.1 


18.9 


0.13 


13 


42 


< 0.8 


< 0.018 








79 


126 


< 6.2 


< 0.035 


1.3 


24.2 


0.89 


84 


51 


< 1.5 


< 0.029 


0.9 


17.0 


0.90 


85 


190 


< 14.0 


< 0.074 


1.0 


17.1 


1.01 


94 


49 


< 0.9 


< 0.019 


0.8 


9.6 


1.44 


125 


98 


< 2.1 


< 0.018 


2.2 


17.5 


2.12 


126 


54 


< 1.3 


< 0.020 


0.8 


6.5 


2.05 


128 


117 


< 1.8 


< 0.018 


1.0 


10.7 


1.63 


130 


27 


< 2.1 


< 0.107 


0.2 


2.3 


1.31 


131 


29 


< 1.4 


< 0.045 




2.6 




133 


128 


< 3.7 


< 0.025 


1.6 


7.2 


3.80 



Note. — Uncertainties iV(N2H+), JV(H2), and JV(CO) are discussed in the text. 
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relative to H2 toward twelve cores, compared with the 
canonical abundance ratio of [C^^O] [H2] = 1.7 x 10^'' 
(jFrerking et al.|[l982D . The CO depletion measurement 
is dependent on the conversion between observed submil- 
limeter continuum emission and H2 column density. As 
described in i i3.2.2| several parameters in this conversion 
are not known to high precision. 

It is also possible that the canonical abundance ratio 
of CO to H2 i nay vary between nearby m olecular clouds. 
For example, iTafalla fc Santiagol (|2004l ) find a best fit 
C^^O abundance relative to H2, X{C^^O) = 0.7 x lO"'^ 
toward the starless core L 1521E in Taurus, which they 
argue is undepleted in CO. lLacv et al.l (|1994D observe di- 
rectly the vibrational absorption of both H2 and CO to- 
ward two infrared sources in Orion, and find [C0]/[H2] = 
2.7 X lO""*, resulting in an abundance ratio for C^*0 
of ~ 4.8 X IQ-^ assu ming [C01/[ C^^01 = 56 ± 25 
(jWilson fc RoodI Il994f ). Using the iLacv et al.l resuh 
then gives /d > 1 for all but three cores, and increases 
the mean C^^O valu e to /d = 8, in agreement with 
IKeto fc Caseliil (I2OO8I) . 

We next compare t he Perseus CO deple tion results 
to the recent work by iChristie et al.l (|2012| ) , who have 
performed the first large-scale, systematic study of the 
CO depletion toward star-forming cores in the Tau- 
rus, Serpens, Ophiuchus, and Orion molecular clouds 
as part of the JCMT Go uld Belt Legacy Survey 
(| Ward-Thompson et al.1 120071 ). They found that while 
cores with high levels of depletion are found within each 
cloud, there is wide variation in the typical CO deple- 
tion found across each cloud. For example, Taurus is 
characterized by relatively high depletion values (mean 
fo = 25 for starless sources), while smaller fd values are 
seen in Ophiuchus (mean fD = 7 for starless sources, and 
fD='i for protostellar sources) , similar to our findings in 
Perseus. The authors suggest the low fu values in Ophi- 
uchus are related to a possible anomalous gr ain size dis- 
tribution in the cloud (iCarrasco et al.l [19731 ). Depletion 
of CO can be suppressed by the presence of small grains, 
because any CO accreted onto the grain will q uickly 
desorb by thermal evaporation (jLeger et al.lll985l ). In- 
triguingly, the presence of anomalous dust emission has 
also been noted in Perseus, where microwave observa- 
tions show evidence of a population of small, spinning 
dust grains in the dense molecular gas on degree scales 
(jWatson et al.l I2005t IPlanck Collaboration et al.l I20H1 ). 
This spinning dust emission is associated, however, with 
the lower density envelopes of the molecular clouds them- 
selves, so it is unclear how the small grains responsible 
for the microwave emission would impact the depletion 
of CO toward the cores themselves. 

Increasing central density sho uld also le ad to 

greater deuteriuni fractionation (iRpueff et all 120051 : 
lAikawa et all 120051 : iCaselli et all [20081). In agreement, 
in tj3.2.1[ we find a significant correlation between in- 
creasing Rr) and central core density for starless cores. 
The deeply embedded sources (those not detected in 
Spitzer catalogs of protostellar sources in Perseus) ap- 
pear to fall along the starless core trend. Protostellar 
cores, in contrast, show a scatter oi Rjj values even at 
high central densities. Comparing Figures [7] and [HI we 
see that the four protostellar sources showing the great- 
est levels of deuterium fractionation have high central 



densities (based on BE modeling), and are also very 
young based on their low Lboi and Ti,oi, while retaining 
significant material in their envelopes (envelope masses 
Menv ^ Mq ) . Figure [H] also shows that no protostellar 
cores have Rjj > 0.05 at i6o;/^i.imm > 3, with only two 
sources showing high Rjy values at Ltoi / Fi.imui > 1. 
Multiple protostellar cores with Lf,oi/^'i.imm < 1 also 
have low Rd, however. The scatter in the Rjj values 
against Lboi and Thoi may be the result of differences 
between individual cores in how the core evolved from 
starless to protostellar. For example, the maximum Rjj 
reached at the starless core stage could vary between 
cores depending on the timescale over which the core 
condensed and collapsed. Nonetheless, the data are con- 
sistent with the general evolutionary scenario where the 
deuterium fractionation in starless cores increases with 
central density, and decreases with time after a proto- 
star has formed and is detected in infrared emission. 

The other parameters are less certain for the Perseus 
cores, and may be the cause of the large scatter in Rd 
in the cor es if they vary acr oss the Perseus molecu- 
lar cloud. ICaselli et al.l ((2008[ ) show that variations in 
the fractional abundance of dust grains will impact the 
deuterium fractionation, where a decrease in the m ini- 
mum radius of dust grains in a iMathis et all ()1977D size 
distribution leads to an decrease in the deuterium frac- 
tionation. We have noted above the possibility that an 
anomalous population of small dust grains in Perseus 
is the cause of the low CO depletion. The deuterium 
fractionation will also decrease with increased cosmic- 
ray ionization rate, as the increased number of electrons 
will dissociate II2D+ before the deuter ium can be trans- 
ferred to other species, such as N2H+ ([Roberts fc Millail 
[2OOOI) . At an angular resolution of several arcminutes, 
the dust temperatures are g reater b y several Ke l vin in 
the east ([Schnee et al.l 120051 ) . which iRidge et al.l ([2OO60 
argue is due to heating from an HIT region driven by 
the nearby early-B star HD 278942 (although note that 
both IC348 and NGC 1333 show localized higher temper- 
atures). On core scales, the gas temperatures determined 
through analysis of the NII3 (1,1) and (2,2) lines are sig- 
nificantly elevated only toward NGC 1333. While the 
elevated dust temperatures on large scales thus do not 
appear to affect the intrinsic core temperatures on small 
scales, they provide evidence of variation in the inter- 
stellar radiation field (ISRF) across the Perseus molecu- 
lar cloud. A significantly stronger ISRF would increase 
the electron abundance, however, and consequently de- 
crease the abundance of molecular ions (such as N2II+) 
through dissociative recombination reactions. As dis- 
cussed in we do not see a systematic variation in 
N2H+ abundances in the Perseus cores across the molec- 
ular cloud. 

Finally, high ratios of ortho-H2 to para-H2 in the gas 
will also lower the abundance of de uterated species by 
decreasing the II2D+ abundance ("Pag ani et al.l [19921: 
Gerlich et al. 2002; Walnislcv et al. 20041: [Flower et al. 
| 200&,Caselh et al.ll2008t iPagani et al.ll2009[: (Pagani et a" 
I2OIII). In contrast to CO, the timescales for some of the 
chemical reactions involved in the deuterium fractiona- 
tion process are longer than the dynamical timescale. In 
particular, H2 forms with an ortho- to para- ratio of 3:1, 
and does not reach the steady-state ratio of 2.7 x 10~^ 
(assuming constant T = lOK and n = lO^'cm"'^) for 



20 



Friesen et al. 



ower "erall[2006h . In regions like Perseus, 
the ortho- to para-H2 ratio, and consequently the deu- 
terium fractionation of species like N2H+, will depend 
on the condensation timescale and thermal history of 
the individual cores. Lower deuterium fractionation val- 
ues may therefore indicate that the cores in Perseus are 
young (t < 10^ yr), in agreement with age estimates by 
Uorgensen et ali ()2008[ ) based on infrared and continuum 
analysis. 



5. SUMMARY 

We have performed a systematic study of the deu- 
terium fractionation of 64 dense cores within the Perseus 
molecular cloud through analysis of N2D"'" observa- 
tions and published N2H+ results, detecting 54 sources 
in N2D+ 2-1, N2D+ 3-2, or both. In the detected 
sources, we find a mean N2H+ deuterium fraction, Rjj = 
A^(N2D+) / 7V(N2H+), of 0.08. We find no significant dif- 
ference in the N2H+ deuterium fraction between starless 
and protostellar cores, and no difference based on core 
location in clustered or isolated environments, although 
isolated cores were least likely to be detected in N2D+ 
emission. 

Kinematically, we find that N2H+ and N2D+ exhibit 
very similar behaviour in the Perseus cores. Both have 
similar, very small, velocity differences with respect to 
NH3 (typically less than ~ 0.07 km s^"'^, and both also 
have similar, slightly larger, velocity differences with re- 
spect to C-'^^O (typically less than ^ 0.19 km s~^); for a 
temperature of 15 K, the sound speed is 0.23 km s~^. The 
small centroid velocity differences between N2D+, N2II+, 
and NH3 suggest that these nitrogen-bearing species all 
occupy a similar volume of each dense core, despite the 
differences in critical densities. 

We compare the deuterium fractionation levels with 
parameters linked to advanced core evolution, and find 
a few significant correlations. First, we find that greater 
N2H+ deuterium fraction correlates with increased II2 
column density and increased central core density, when 
considering the entire core sample. The trend with 
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A^(Il2) is stronger when limited to only protostellar cores, 
while the trend with log(nc) is stronger when limited to 
only starless cores. We additionally find that the N2H"'" 
deuterium fraction decreases with increased Tboi for pro- 
tostellar sources. Significant correlations were not found 
between Rd and other parameters, including N2II+ line 
width, kinetic gas temperature, or protostellar bolomet- 
ric luminosity. The data are nevertheless in agreement 
with the general expectation that the N2II+ deuterium 
fraction will decrease toward protostellar cores as the em- 
bedded protostar warms its immediate environment 

We show that the Perseus cores are characterized by 
low CO depletion values relative to previous studies of 
star forming cores, with a mean fu of 3.1. The CO deple- 
tion values are similar, however, to those found for cores 
in the Ophiuchus molecular cloud. We suggest that the 
low average CO depletion factor is the dominant mech- 
anism that constrains the average deuterium fractiona- 
tion in the Perseus cores to small values. While current 
equilibrium and dynamic chemical models are able to re- 
produce the range of deuterium fractionation values we 
find in Perseus, reproducing the scatter across the cores 
requires variation in parameters such as the ionization 
fraction or the ortho- to para-H2 ratio across the cloud, 
or a range in core evolution timescales. 
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